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I .  INTRODUCTION 


1  2 

In  recent  work,  one-  and  two-dimensional  calculations  of  the  properties 
of  arcs  in  arc-driven  rail  guns  were  carried  out.  The  principal  purpose  of 
that  work  was  two-fold.  First,  we  wanted  to  determine  the  most  reasonable  arc 
materials  and  whether,  by  a  judicious  choice  of  material,  we  could  control  the 
arc's  thermodynamic  properties.  Second,  we  wanted  to  provide  the  basis  for  in¬ 
vestigating  what  erosive  effects  the  arc  had  on  the  rails  and  projectile.  In 
the  present  work  we  undertake  a  preliminary  investigation  of  what  the  erosive 
effects  are. 


Specifically,  we  wish  to  provide  a  simple  model  for  determining  thermal- 
energy  generation  in  and  transfer  to  the  rails,  and  for  assessing  the  damage 
to  the  rails  caused  by  that  thermal  energy.  The  energy  generated  in  the  rails 
will  be  assumed  to  arise  from  ohmic  heating,  while  that  transferred  from  the 
arc  to  the  rails  will  be  assumed  to  be  radiative.  Although  extensive  experi¬ 
mental  work  designed  to  examine  carefully  rail  damage  has  not  been  carried  out, 
early  work  has  indicated  the  greatest  damage  near  the  breech  end  of  the  gun, 
but  practically  none  farther  down  the  gun  tube.  It  is  not  clear,  however,  that 
even  at  the  breech  end,  the  damage  is  thermal  in  origin  and  it  is  important  to 
determine  theoretically  under  what  conditions  such  damage  might  be  significant. 


The  model  used  in  Refs.  1  and  2  is  shown  schematically  in  Fig.  1.  Ini¬ 
tially,  the  rails  (the  upper  and  lower  sides  in  the  figure),  which  are  infi¬ 
nitely  extended  in  the  z  direction,  are  supplied  a  pulsed  current  i^  which  is 

assumed  to  be  constant  in  time  and  which  flows  in  the  direction  indicated.  In 
our  previous  work  the  rails  were  assumed  to  be  perfect  conductors  so  that  the 
current  was  carried  only  along  their  inner  surface  and  their  thickness  was  in¬ 
significant;  in  the  figure,  however,  we  have  indicated  the  finite  thickness  d 
of  the  rails.  For  finite  conductivity,  of  course,  the  current  is  carried  with¬ 
in  some  skin  depth  A  from  the  rail  surface.  The  current  is  also  conducted 
through  the  arc  which  has  length  i  .  The  resulting  current  configuration  then 

produces  a  magnetic  field  between  the  rails  in  the  space  bounded  by  the  power 
supply  and  the  arc,  as  well  as  in  the  arc  itself.  This  induction  field,  in 
turn,  interacts  with  the  current  through  the  arc  producing  a  high-pressure 
plasma  which  accelerates  the  projectile,  shown  by  the  shaded  area,  in  the  x  di¬ 
rection.  One-  and  two-dimensional  calculations  were  undertaken  to  determine 
the  properties  of  the  arc  under  the  assumption  that  those  properties  were  steady 
as  seen  by  an  observer  in  the  accelerating  reference  frame.  Some  justification 
for  the  steady-state  approximation  was  given  for  the  acceleration  times  of  in¬ 
terest  . 


In  the  present  calculation  we  wish  to  consider  the  problem  of  thermal - 
energy  transfer  from  the  arc  to  the  rails.  For  completeness,  we  will  also  ac¬ 
count  for  resistive  heating  in  the  rails,  though  in  most  practical  problems 


2J.D.  Powell  and  J.H.  Batteh ,  '"Plasma  Dynccnies  of  an  Arc-Driven,  Electromag¬ 
netic,  Projectile  Accelerator J .  Appl.  Phys.  52,  2717  (1981).  See  also 
"Plasma  Dynamics  of  an  Arc -Driven  Rail  Guny  rr  Ba  tlistic  Research  laboratory 
Report  No.  ARBRL-TR-02267 ,  September  1980  (AD  A092345). 

D  Powell  and  J.H.  Batteh,  "Two-Dimensional  Plasma  Model  for  the  Arc-Driven 
Rail  Gun,"  J.  Appl.  Phys.  54,  2242  (1983).  See  also  J.D.  Powell,  "Two- 
Dimensional  Model  for  Arc  Dynamics  in  the  Rail  Gun, "  Ballzstzc  Research  labora¬ 
tory  Report  No.  ARBRL-TR-02423,  October  1982  (AD  A120046). 
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Figure  1.  Rail  gun  model. 


resistive  heating  is  insignificant  compared  to  heating  from  the  arc,  especially 
if  the  rail  thickness  is  comparable  to  or  larger  than  the  electrical  skin  depth. 

Consider  the  point  in  Fig.  1.  When  the  leading  edge  of  the  arc  reaches 

the  point,  the  temperature  there  will  rise  both  because  of  heat  transfer  from 
arc  to  rail,  as  well  as  because  of  ohmic  heating  initiated  in  the  rail  once  cur¬ 
rent  begins  to  flow.  Once  the  arc  has  completely  passed  x^,  radiative  transfer 

from  arc  to  rail  is  assumed  to  be  small,  though  the  temperature  in  the  interior 
of  the  rail  still  rises  because  of  ohmic  heating.  If  and  when  the  temperature 
at  the  rail  surface  reaches  the  melting  temperature,  T  ,  the  temperature  there 

will  remain  constant  but  the  surface  will  recede  from  its  initial  position  at 
y  =  0.  Therefore,  the  location  of  the  rail  surface  is  a  time-dependent  quantity 
given  by  y  =  s(t)  ,  and  the  magnitude  of  s(t)  is  indicative  of  the  mass  lost  at 
the  gun-tube  surface  due  to  melting.  We  want  to  calculate  the  location  of  the 
surface  as  well  as  the  temperature,  magnetic  induction,  and  current  density  with¬ 
in  the  rail . 

The  calculation  is  characterized  by  a  number  of  assumptions  made  to  make 
the  analysis  as  simple  as  possible.  Perhaps  the  most  significant  is  that  the 
arc-dynamic  problem  can  be  uncoupled  from  the  rail-heating  problem.  It  is, 
for  example,  assumed  that  the  radiation  flux  from  the  arc  incident  on  the  sur¬ 
face  of  the  rails  is  independent  of  time  and  unaffected  by  the  rise  in  the 

rail  temperature.  As  has  been  argued  before,*  such  an  assumption  is  clearly 
reasonable  since  the  arc  temperature  is  more  than  an  order  of  magnitude  higher 
than  that  of  the  rails.  The  arc  serves,  therefore,  essentially  as  a  high- 
temperature  reservoir  whose  thermodynamic  properties  are  nearly  steady,  while 
those  of  the  rails  are  quite  unsteady.  It  is  also  assumed  that  the  magnetic 
flux  density  at  the  rail -arc  interface  varies  linearly  from  zero  at  the  front  of 
the  arc  to  yj  at  the  back  of  the  arc.  Here  y  is  the  magnetic  permeability  and 
j  is  given  by 

d 

j  =  /  J*(xo *>')  dy*  n.n 

0 

where  J  is  the  current  density  in  the  rails  and  d  is  the  rail  thickness.  That 
the  induction  field  should  take  on  these  values  in  front  of  and  behind  the  arc 
can  be  proved  easily  by  Ampere's  law.  While  the  linear  variation  at  the  in¬ 
terface  is  not  strictly  correct  for  the  problem  at  hand,  it  would  be  valid  if 
the  rails  were  perfect  conductors  and  if  the  conductivity  of  the  arc  were  con¬ 
stant,  and  it  is  nearly  valid  for  more  realistic  values  of  the  arc  conductiv- 
1  2 

ity.  ’  For  rails  which  are  not  perfect  conductors,  determining  the  value  of 
the  induction  field  at  the  interface  would  necessitate  solving  the  magnetic 
diffusion  equation  in  both  the  rails  and  arc,  and  applying  the  appropriate 
boundary  conditions  at  the  interface.  For  such  a  calculation  the  electromag¬ 
netic  properties  of  the  arc  are  coupled  to  those  of  the  rails.  We  expect, 
however,  that  since  the  rail  conductivity  is  several  orders  of  magnitude 
higher  than  that  of  the  arc,  the  assumption  of  perfect  conductivity  is  satis¬ 
factory  for  determining  the  interface  boundary  condition.  Furthermore,  any 
errors  made  in  such  an  assumption  should  be  insignificant  since,  in  practical 
problems,  while  the  arc  is  in  contact  with  the  point  x0>  the  dominant  mechanism 

for  heat  production  there  is  radiation  from  the  arc. 
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The  second  principal  assumption  is  that  all  the  radiative  flux  from  the 
arc  that  is  incident  on  x^  is  absorbed  by  the  rails,  just  as  if  the  arc  were 

radiating  into  a  vacuum.  No  doubt  some  of  this  incident  energy  is  reflected 
back  into  the  arc,  but  it  is  difficult  to  estimate  how  much.  Furthermore,  by 
assuming  total  absorption,  we  should  overestimate  the  erosive  effects  on  the 
rails  and  can  easily  alter  the  heat  flux  when  more  reliable  estimates  arc  known. 
The  value  of  the  heat  flux  incident  on  the  rails  is  assumed  to  be  given  by 

Q=2°ST*  (1.2) 

where  T,  is  the  arc  temperature  at  the  boundary  and  where  o  is  Stefan’s  con- 

°  -  8  J  ** 

stant,  namely,  5.67  x  10  ^ -  .  This  approximate  expression  has  been 

j  2  m  K  s 
discussed  elsewhere.  ’ 

The  third  major  assumption  is  that  any  melted  portion  of  the  rails  is 
"swept  away"  by  the  motion  of  the  arc  so  that  the  radiation  flux  is  always  in¬ 
cident  on  the  solid  rail  surface.  It  might  appear  that  such  an  assumption 
would  be  unreasonable  because  one  would  expect  that  the  liquid  rail  material 
would  remain  on  the  surface  because  of  inertial  effects  during  the  time  of  in¬ 
terest  (tens  of  microseconds).  However,  for  the  problem  at  hand,  the  pressures 
at  the  rail  surface  are  of  the  order  of  hundreds  of  atmospheres  and  these  large 
pressures  may  indeed  produce  large-scale  motion  during  the  time  of  contact. 

At  any  rate,  the  assumption  should  again  overestimate  the  erosive  effects  since 
none  of  the  energy  is  expended  in  heating  the  liquid  or  gaseous  phase. 

Three  additional  assumptions  in  the  calculation  merit  brief  mention. 

First,  we  will  neglect  diffusion  within  the  rail  of  both  the  magnetic  induc¬ 
tion  field  and  the  temperature  in  the  direction  of  propagation,  and  assume 
that  it  occurs  only  in  the  y  direction.  The  assumption  is  likely  to  be  valid 
for  projectile  velocities  that  are  high  compared  to  the  speed  with  which  both 
quantities  can  diffuse  into  the  rail.  Such  is  the  case  in  any  practical  ex¬ 
periment.  Second,  we  will  assume  that  the  major  form  of  energy  transport 
within  the  rails  is  ordinary  heat  conduction  and  will  neglect  radiation.  That 
radiation  is  negligible  below  the  melting  temperature  of  rail  materials  is 
easily  demonstrated.  Finally,  we  neglect  the  temperature  dependence  of  the 
electrical  and  thermal  conductivities  and  of  the  specific  heat  of  the  rails. 
This  assumption  is  made  partly  because  reliable  expressions  are  not  known 
everywhere  in  the  temperature  range  under  consideration,  and  partly  because 
analytic  solutions  can  be  obtained  in  special  cases  and  these  provide  a  useful 
comparison  with  the  numerical  solutions. 


II.  FORMALISM 

Consider  point  x^  in  Fig.  1  and  note  that  prior  to  acceleration  the  lead¬ 
ing  edge  of  the  arc  is  located  at  x  =  0.  Then  since  the  acceleration  a  is  con¬ 
stant,  the  arc  velocity  when  the  leading  edge  reaches  x^  is  given  by 
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vo  =  ^ 2a  xo3 


(2.1) 


1/2 


If  we  denote  this  time  by  t  =0,  we  can  write  that  in  time  t  the  leading  edge 
will  move  beyond  xQ  by  a  distance  D  given  by 


D  =  v^t  +  1/2  a.  t 


2 


(2.2) 


Furthermore,  since  it  is  assumed  that  the  magnetic  induction  field  varies 
linearly  from  zero  to  pj  as  one  moves  from  the  front  to  the  rear  of  the  arc, 
we  clearly  have  at 


B(s,t) 


+  1/2  a  t2) 


‘1*0 


(2.3) 


B (s  ,t)  =  pj 


t  >  t 


0 


Here  8.^  represents  the  total  length  of  the  arc  and  tQ  is  the  time  for  which 
the  arc  is  in  contact  with  the  point  x^.  One  can  easily  show  that 


t0  =  /27*  (/xo  +  -  ’'*0  3 


(2.4) 


Given  the  time -dependent  boundary  condition  at  y  =  s,  the  B  field  for  y  >  s 
can  be  determined  from  the  magnetic  diffusion  equation 


pa 


3B  32B 


(2.5) 


where  diffusion  in  x  direction  has  been  neglected  to  accord  with  the  assump¬ 
tions  in  Sec.  I.  The  current  density  then  follows  from  the  relation 


J 


1_  9B_ 

‘  J  X  ~  p  3y  * 


(2.6) 


Heat  conduction  in  melting  solids  has  been  considered  by  various  authors. 
Probably  the  treatment  most  closely  related  to  the  present  work  is  that  by 

Landau-5  to  which  the  reader  is  referred  for  additional  references.  Generally 
we  can  write  the  energy-conservation  relation  within  the  rail  as 


p  c  It  +  iy  =  j2/a  ■  EL 


(2.7) 


7 

H.G.  Landau,  "Heat  Conduction  in  a  Melting  Solid,"  Quart.  Appl.  Math.  8^  81 
(1950). 


11 


where  p  represents  the  solid-phase  density,  C  is  the  specific  heat,  T  the  tem 
perature,  q  the  heat  flux,  o  the  electrical  conductivity,  and  E^  the  energy 

absorbed  per  unit  volume  and  time  by  the  solid  in  the  melting  process.  If  we 

2 

identify  the  term  J  /a  as  the  energy  dissipated  per  unit  volume  and  time  aris 
ing  from  resistive  heating,  the  physical  meaning  of  Eq.  (2.7)  is  apparent. 

Now  if  the  solid-liquid  transition  takes  place  at  a  well-defined  tempera 

ture  T  as  for  most  metals,  we  can  write 
m 

El  =  Lp  «<T-Tb)  |1  (2.8) 

where  6  is  a  Dirac  delta  function  and  L  is  the  heat  of  fusion,  i.e.,  the  heat 
absorbed  per  unit  mass  of  solid  to  produce  an  equivalent  amount  of  liquid. 
Furthermore,  from  properties  of  the  delta  function, 

■5(T'T")  ’  I  St  ,  tty  - s)  •  t2-9) 

'  3y  '  v=s 

It  is  evident  on  physical  grounds,  however,  that 

(§p>  <  0  (2.10) 

y-s 

and  that  at  y  =  s,  T  changes  only  by  virtue  of  the  moving  boundary.  Conse¬ 
quently,  it  follows  that 


( —  ds  ,-3T ^ 

l3tJ  "  '  dt  l3yJ 
y=s  y=s 

and  Eq.  (2.7)  finally  becomes 

.2 


PC 


3T  _  J_ 
3t  +  3y  ”  a 


lp  «(y 


-  s) 


(2.11) 


(2.12) 


The  last  term  in  Eq.  (2.12)  is  zero  except  at  y  =  s  and  the  boundary 
conditions  to  which  it  gives  rise  can  be  obtained  by  integrating  Eq.  (2.12) 
across  a  negligibly  small  interval  centered  at  y  =  a  •  We  find 

q  (s  +  e)  -  q  (s  -  £)  =  -  Lp  .  (2.13) 

At  s  -  c  the  flux  is  just  that  incident  from  the  arc,  namely, 
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Q  =  2  O  T?  ; 
<  s  b  ’ 


(2.14) 


at  s  +  e  we  have,  assuming  that  ordinary  heat  conduction  is  dominant. 


q  =  -  k 


3T 

97 


(2.15) 


where  k  is  the  rail  thermal  conductivity.  Thus,  Eq.  (2.13)  becomes 


k  ft 

l3yJ 

y=s 


Lp  57 


(2.16) 


Equations  (2.5),  (2.6),  (2.12)  and  (2.16)  are  now  sufficient,  with  the 
boundary  conditions  discussed,  to  determine  the  temperature,  magnetic  induc¬ 
tion,  current  density,  and  boundary  location.  It  must  be  remembered,  however, 
that  melting  will  temporarily  cease  if  T(s,t)  drops  below  T  ,  and  in  that 


ds 


It  is  convenient  to  transform  the  governing  equations 
to  a  frame  of  reference  that  moves  with  the  melting  boundary,  so  we  let 


case,  becomes  zero 


n  =  y  -  s(t) 


(2.17) 


and  obtain 


3B  32B  ds  3B 

y0  5t  =  +  ya  dt  37 

.  1  3B 

J  ‘  u  3n 


„  3T  J  3  ,,  3T.  _  ds  3T 

pc  3t  =  -  +  37  (k  an5  +  pc  - 


9n 


dt  3n 


Q  -  -  k  ^n=0  +  3F 


(2.18) 

(2.19) 

(2.20) 
(2.21) 


These  equations  must  be  solved  subject  to  the  boundary  and  initial  conditions: 
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8(0,0  =  p-  CvQ  +  1/2  a  t2)  t  <  tQ 

3. 

8(0, t)  =  yj  t  >  tQ 

B(n,0)  =  0  (2.22) 

B (d-s , t)  =  0 

$W.  * 0 

T(n.O)  -  TQ 

3F>0“T(0-t)  -T. 

where  Tq  is  the  ambient  temperature. 


III.  CALCULATIONS 
A.  Special -Case  Analytic  Solutions 


We  will  be  concerned  primarily  in  this  section  with  obtaining  numerical 
solutions  to  Eqs.  (2 . 18) -(2 . 21)  .  Prior  to  undertaking  the  numerical  calcula¬ 
tions,  however,  it  is  useful  to  carry  out  some  special-case  analytic  solutions 
for  checking  the  computer  program  as  well  as  giving  some  physical  insight  into 
the  nature  of  the  solution  when  special  conditions  apply.  It  will  be  assumed  in 
this  section  that  the  rails  are  infinitely  thick,  an  assumption  that  is  reason¬ 
able  provided  the  rail  thickness  d  is  large  compared  to  the  electrical  skin 
depth.  This  condition  is  satisfied  in  all  cases  that  will  be  considered. 

Preparatory  to  carrying  out  the  analytic  solutions,  it  is  useful  to  point 

4 

out  the  following  general  results.  Consider  the  one-dimensional  diffusion 
equation  of  the  form 

-  a2  ||  =  -  4ir  F  (n,  t)  (3.1) 

3n 


4 

P.M.  Morse  and  H.  Feshbaoh,  Methods  of  Theoretical  Physics  (McGraw-Hill,  New 
York,  1953),  Chap.  7. 
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which  obeys  appropriate  boundary  and  initial  conditions  and  where  is  defined 
in  the  half -space  r)  _>  0 .  The  general  solution  of  Eq.  (3.1)  may  be  written 


t  00  2  00 

<Kn,t)  -  /  dt'  /  drr  F(rT,  t')  G(n,tjn',0  +  ^  /Q  diT  Wt^=0 


0  0 
t 


-W  '  dt'  'G  ST  -  *  Iff-  v.  0 


(3.2) 


where  G  is  the  appropriate  Green’ s  function.  For  a  semi -infinite  medium  the 
Green’s  function  can  be  constructed  from  those  appropriate  for  an  infinite 
medium.  One  finds 


G(n,t;rT,0 


exp  [ 


-a2(n-n") 

~ *Tt-t') 


exp 


r-ct  (n+n ") 

l— Xt-t  O' 


2 

-1} 


(3.3) 


3C 

where  the  upper  sign  is  to  be  taken  if  is  to  vanish  at  ri'  =  0  (homogeneous 

Neumann  condition)  and  the  lower  sign  if  G  is  to  vanish  at  r\'  =  0  (homogene¬ 
ous  Dirichlet  condition) .  These  results  will  prove  useful  in  obtaining  the 

special -case  solutions. 


1 .  Current  Density  and  Magnetic  Induction  Prior  to  Melting 


Analytic  expressions  for  both  the  current  density  and  the  induction 
field  can  be  obtained  prior  to  the  onset  of  melting  from  Eqs.  (2.18)  and  (2.19) 
and  from  the  boundary  conditions  in  Eqs.  (2.22).  The  last  term  in  Eq.  (2.18) 
is,  of  course,  set  equal  to  zero.  Using  Eqs.  (3.1)-(3.3)  and  restricting  our¬ 
selves  to  times  less  than  t^,  we  obtain 


B(n.t)  = 


1/2  t 

HHL  (H2.)  f 

2i  1  J 


(vo 


u 


t'+  1/2  a  t 

:t  -  to3'* 


.2 


r  V  <7  n  1  J  + 

exP  4(t-tJ)]  dt 


(3.4) 


The  integral  in  Eq.  (3.4)  can  be  evaluated  by  making  the  change  of  variable 


2 

r 


2 

pan 

4(t  -  t') 


(3.5) 


and  using  the  easily  derived  reduction  formula 


-r 


2j 


j _  ~  e  2 _ [ 

dr  =  2j -1  (2j  -T)  J 


-r 


r  -  (2j  -  1)  r 

One  finds  after  considerable  algebra 


2j  “  2 


dr 


(3.6) 
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B(n,  t)  =  Hi  (v0t  +  1/2  a  t2)  Erfc  (r0) 
a 


2  2  2  3  ^  * 

,  2  y  ajq  s  a  pari  _-rj;  .  v0  +  at,  re"r0 

1  x  '  IS"  4  ^  1— 

a  0 

-  /n  Erfc  (rn)]} 


(3.7) 


where 


2  _  yon' 


Tn  = 


4t 


(3.3) 


and  where  Erfc  is  the  complementary  error  function 


Erfc  (r  )  =  —  /  e-r  dr  . 
°  ft  r„ 


(3.9) 


The  current  density  then  follows  by  direct  differentiation  in  accordance  with 
Eq .  (2 . 19) .  We  find 


1/2 


J(n,t)  =  -  {  (vn  +  at)[(T^)  e~r0  ~Sn  q  Erfc  (rn)] 

l  Sn 
a 


0  — '  l'yo 

2 


ayo  e"r0  2q2  f4t^^2  -r2  +  2/ffq^  r  ^  f 

—  “3 - 3“W  6  °  - S-  Erfc  (rQ)] } 


(3.10) 


2 .  Temperature  Profile  Prior  to  Melting  for  Negligible  Ohmic  Heating 

As  indicated  previously,  ohmic  heating  is  frequently  negligible  in 
cases  of  practical  interest.  In  such  a  case,  Eq.  (2.20)  can  be  solved  prior 
to  the  time  of  melting  via  the  techniques  discussed  previously.  Equation 

8T 

(2.21)  is  used  to  supply  the  boundary  condition  on  5 —  at  q  =  0  and  the  "raelt- 

dq 

ing  term"  in  both  Eqs.  (2.20)  and  (2.21)  is  set  equal  to  zero.  We  find  for 
times  t  <  t,_ 

2 

pCqZ 

(3.11) 


0 


TCn.O  =t 0+-a-  J* 

0  /froCF  0  (t  -  t')1/2 


The  integral  can  be  evaluated  as  before  by  using  Eqs.  (3.5)  and  (3.6)  and  we 
find 


1/2  2 

T(n,t)  =t0.3[§  exp  (  )  -  n  Erfc  n)] 


(3.12) 
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From  Eq.  (3.12)  an  approximate  time  to  melting  t  can  be  derived  by  setting 

ri  =  0  and  T(n,t)  =  T  .  We  obtain 

m 


t 


m 


i  p  C  k  (T 
K  v  m 


(3.13) 


3.  Steady-State  Melting  for  Negligible  Ohmic  Heating 

For  sufficiently  long  times,  one  might  expect  that  Eqs.  (2.18) -(2.21) 
would  have  a  steady-state  solution  provided  the  resistive-heating  term  can  be 
neglected.  Assuming  such  a  solution,  we  obtain  from  Eqs.  (2.18)  and  (2.20) 


d2B  dB  _ 

— =-  +  uo  v  j—  =  0 
,2  H  m  dn 
dn 


and 


.  d2T  „  dT 

k  — T  +  P  C  V  -r— 

.2  in  dn 

dn 


(3.14) 


(3.15) 


where 


(3.16) 


is  the  (constant)  propagation  velocity  of  the  melting  surface. 
Equations  (3.14)  and  (3.15)  have  the  obvious  solutions 


T  -  <Tm  -  V  e'pC  vm  n/k  ♦  Tq 

Furthermore,  from  (2.21)  and  (3.18) 
velocity 


(3.17) 

(3.18) 

we  obtain  for  the  steady-state  melting 


Q  .  (3.19) 

Vm  '  PC  (Tm  -  Tq)  ”  Lp 

The  nontrivial  solution  arises  when  Q,  the  constant  flux  radiated  from  the  arc, 
is  nonzero.  Therefore,  one  can  expect  such  a  steady  solution  to  be  achieved 
provided  that  the  arc  is  in  contact  with  xQ  long  enough  for  a  steady  state  to 
be  achieved. 
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B.  Numerical  Analysis 


To  solve  Eqs.  (2 .18) -(2.21) ,  we  divided  the  plane  into  a  grid,  re¬ 

presented  the  derivatives  by  second-order  finite  differences,  and  solved  by  the 

standard  implicit  technique.  Not  only  were  numerical  calculations  carried 
out  for  cases  of  general  interest,  but  also  for  cases  corresponding  to  the 
special-case  analytic  solutions  in  Sec.  Ill  A.  We  could  thereby  check  the 
validity  of  the  computer  program  as  well  as  determine  suitable  step  sizes  Ari 
and  At.  For  all  cases  undertaken,  we  found  acceptable  values  of  Aq  and  At  of 

-6  $ 

the  order  of  5  X  10  m  and  5  X  10  s,  respectively,  while  the  arc  was  in  con¬ 
tact  with  the  point  x^;  some  time  after  the  arc  had  passed,  however,  and  radia¬ 
tion  was  assumed  negligible,  we  were  able  to  increase  At  by  an  order  of  magni¬ 
tude  or  more. 

We  have  performed  calculations  in  which  we  both  neglected  and  included 
the  resistive-heating  losses.  We  consistently  found,  however,  that  for  all 
cases  of  interest  to  us  (very  high-temperature  arcs  and  rails  whose  thickness 
was  greater  than  or  comparable  to  a  skin  depth) ,  the  resistive  contributions 
were  utterly  negligible.  Therefore,  in  the  discussion  which  follows,  the  re¬ 
sistive  terms  have  been  neglected.  Doing  so  immensely  facilitates  numerical 
calculations  because  the  speed  at  which  the  magnetic  field  diffuses  into  the 
rail  is  much  larger  than  the  speed  at  which  the  temperature  diffuses.  Conse¬ 
quently,  when  considering  the  negligible  resistive  losses,  it  was  necessary  to 
employ  far  more  grid  points  in  the  calculation  than  was  necessary  in  the 
simpler  case. 

C.  Analysis  of  Jamison-Burden  Experimental  Results 

In  a  recent  experiment,  Jamison  and  Burden  (JB)  have  employed  a  rail  gun^ 
made  of  copper  rails  with  the  projectile  driven  by  a  copper-vapor  arc  to  ac¬ 
celerate  a  2.5g  projectile  to  about  a  kilometer  per  second.  Values  of  the  cur¬ 
rent,  breech  voltage,  velocity,  acceleration,  and  arc  length  were  measured  at 
several  points  along  the  gun  tube.  These  data  were  then  used  as  input  para¬ 
meters  in  the  numerical  solution  of  Eqs.  (2.20)  and  (2.21),  [The  resistive 
contribution  to  the  temperature  rise  was  found  to  be  negligible  in  this  experi¬ 
ment.  Therefore,  J  was  set  equal  to  zero  in  Eq.  (2.20)  and  Eqs.  (2.18)  and 
(2.19)  were  not  solved.]  In  doing  such  an  analysis  of  course  it  is  necessary 
to  assume  that  all  significant  parameters  associated  with  the  rail  gun  change 
in  a  quasi -static  manner  since  the  original  derivations  were  for  steady-state 
operation. 

In  Table  I  are  shown  the  material  properties  of  copper  used  in  the  calcu¬ 
lation  as  well  as  the  experimental  data  taken  at  two  points  along  the  gun  tube, 
at  xQ  =  0.156m  and  xQ  =  0.51m.  The  first  five  entries  in  the  columns  headed 


^ B .  Carnahan,  H.A.  Luther,  and  J.O.  Wilkes,  Agglied  Numerical  Methods.  (Wiley, 
New  York,  1969),  Chap.  7. 

£ 

K.A.  Jamison  and  H.S.  Burden,  "A  Laboratory  Arc-Driven  Rail  Gun,"  Ballistic 
Research  Laboratory  Report  Noi  ARBRL-TR-Q2502,  June  1982  (AD  A1S11SS). 
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TABLE  I. 

Experimental  Data  and  Input 
Calculations 

Parameters  for  Numerical 

Material  Propert ies 

of  Rails 

Quant ity 

Value 

k 

400  watts  /(m  K) 

0 

5.8  x  10  mho/m 

L 

2.1  x  io5  .T/kg 

C 

390  .T/(K  kg) 

P 

8.9  x  io3  kg/m3 

T 

m 

Data 

1360  K 

Quant ity 

Values  at  x^  =  0.156m 

Values  at  x^  =  0.51m 

vo 

S60  m/s 

757  m/s 

a 

5.8  x  io5  m/s2 

1.69  x  ioS  m/s2 

j 

3.94  x  io6  A/m 

1.97  x  io6  A/m 

vo 

156  volts 

132  volts 

a 

a 

3 . 3  cm 

6.8  cm 

*0 

58  ps 

90  }is 

Tb 

1 .69  x  104  K 

1.14  x  104  K 

"DATA"  were  taken  directly  from  the  experimental  results7  and  arc  self- 
explanatory.  It  should  be  mentioned,  however,  that  the  muzzle  voltage  V^  was 

not  directly  measured  in  the  experiment.  The  breech  voltage  was  measured,  how¬ 
ever,  and  the  muzzle  voltage  estimated  by  measuring  and  subtracting  the  induc¬ 
tive  contribution  to  the  breech  voltage.  The  last  two  entries  in  the  table, 
needed  in  the  numerical  calculations,  were  calculated  for  each  point  in  terms 
of  the  experimental  data.  The  parameter  tR  was  calculated  by  using  Cq .  (2.4), 

with  the  actual  initial  velocity  v^  substituted  for  the  quantity  ►Tox^  ,  and 

the  boundary  temperature  T^  from  the  approximate  expression  derived  in  the 
Appendix: 


1/4 


^b  ^4a  l 


(3.20) 


Jamison  and  II. S. 


Burden 


(private  communication) . 
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Results  of  the  numerical  calculation  at  the  point  xQ  =  0.156m  are  shown 

in  Figs.  2-5.  In  Fig.  2  are  four  temperature  profiles.  Profile  A  shows  the 
temperature  distribution  as  a  function  of  depth  just  prior  to  the  time  melting 
begins.  The  temperature  at  n  =  0  is  about  1335  K  and  the  arc  has  been  in  con¬ 
tact  with  Xq  for  12  ys .  A  nonnegligible  rise  in  the  rail  temperature  can  be 

seen  for  about  100  pm  beyond  the  rail  surface  at  n  =  0.  Profile  B  shows  the 
temperature  distribution  within  the  rail  at  time  t  =  57  ys.  As  can  be  seen 
from  the  value  of  tp  in  Table  I,  this  time  is  just  prior  to  the  time  that  the 

arc  leaves  x_.  The  temperature  at  p  =  0  is,  of  course,  T  =  1360  K.  The  ex- 
u  in 

tent  to  which  Profile  B  is  approximated  by  the  steady-state  result  can  be  seen 
by  comparing  Profiles  B  and  C,  since  Profile  C  is  just  a  graph  of  Eq.  (3.18). 

As  can  be  seen,  the  steady-state  result  has  nearly  been  reached.  Finally,  Pro¬ 
file  D  shows  the  temperature  distribution  at  t  =  90  ys,  i.e.,  long  after  the 
arc  has  passed.  Melting  has  ceased  as  would  be  expected  and  the  temperature 
at  the  surface  has  fallen  by  nearly  a  factor  of  two.  It  should  be  noted  that 
in  terms  of  the  original  location  of  the  rail  surface,  curves  B  and  D  have 
been  shifted  by  about  60  ym  (see  Fig.  5)  because  of  rail  melting. 

In  Fig.  3  is  plotted  the  temperature  of  the  rail  surface  (located  at  p  = 

0)  as  a  function  of  time.  We  observe  a  steady,  rapid  rise  in  the  temperature 
for  about  the  first  fourteen  microseconds  until  the  melting  temperature  Tm  is 

reached.  The  time  to  melting  can  be  predicted  analytically  from  Eq.  (3.13). 
Subsequently,  the  temperature  at  the  surface  remains  constant  until  the  arc 
has  passed  (t  'V  59  ys) .  Continuous  cooling  then  follows  as  heat  is  conducted 
into  the  bulk  of  the  rail 

In  Fig.  4  is  shown  the  velocity  of  recession  of  the  melting  surface,  v  , 

as  a  function  of  time.  We  see  a  rapid  rise  in  the  velocity  during  the  early 
stages  of  melting  followed  by  a  slower  approach  to  the  steady  value  of  1.7 
m/s,  predicted  by  Eq.  (3.19).  Just  prior  to  the  time  the  arc  has  left  xn  and 

melting  has  ceased,  the  surface  velocity  has  approached  about  85%  of  the  steady- 
state  value . 

Finally,  in  Fig.  5,  we  show  the  position  of  the  rail  surface  relative  to 
its  initial  position  (the  parameter  s  defined  in  Sec.  I)  as  a  function  of  time. 
As  can  be  seen,  between  fifty  and  sixty  microns  of  the  rails  will  have  melted 
during  the  time  in  which  the  arc  was  in  contact  with  x^.  After  the  arc  has 
passed  no  further  melting  occurs. 

At  the  point  xQ  =  0.51m,  to  which  the  second  set  of  data  in  Table  I  ap¬ 
plies,  we  have  from  Eq.  (2.4),  t^  =  90  ys .  From  Eq.  (3.13),  however,  the  time 
to  melting  t  is  332  ys  so  we  conclude  that  no  melting  of  the  rails  occurs  at 

-  .  .  7*  m 

this  point . 

D.  Other  Calculations 

The  JB  experiment  is  the  only  case  of  which  we  have  sufficient  experimen¬ 
tal  data  to  allow  us  to  calculate  rail-melting  phenomena  using  only  measured 
quantities.  To  get  some  idea  of  what  the  effects  might  be  in  other  guns,  how¬ 
ever,  we  have  used  results  calculated  in  our  previous  work  as  input  for  our 
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DISTANCE 

Figure  2.  Temperature  profiles  as  a  function  of  distance  behind  the  melting  surface  (n) ;  A,  at  t  =  12  ps 
(just  prior  to  melting);  B,  at  t  =  57  ps  (just  prior  to  the  cessation  of  melting);  C,  steady- 
state  melting  profile;  D  .at  t  =  90  ps  (long  after  the  arc  has  passed).  The  distance  scale  is 
in  microns . 
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Figure  3.  Surface  temperature  vs.  time. 


The  time  scale  is  in  microseconds. 
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Figure  5.  Position  of  solid  surface  relative  to  initial  position,  s, 
in  microns,  vs,  time,  in  microseconds. 
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thermal -energy  transport  model.  Surely  the  greatest  source  of  uncertainty  in 
these  data  is  for  the  arc  length  SL^.  Indeed,  this  quantity  was  unknown  in  our 

previous  work  and  a  value  of  the  order  of  10  cm  was  assumed,  as  had  been  done 

g 

in  previous  arc-dynamic  calculations. 

Input  data  and  results  for  two  calculations  are  shown  in  Table  II.  'the 
first  is  for  a  gun  comparable  in  size  to  that  used  by  Rashleigh  and  Marshall 

9 

(RM)  in  their  well-known  experiment,  and  the  second  was  for  a  larger  gun  com¬ 
parable  to  that  recently  constructed  by  Deis,  McNab,  and  coworkers.10  Two- 

2 

dimensional  calculations  to  determine  properties  of  the  arc  in  both  these  guns 
were  undertaken  in  our  previous  work. 


TABLE  II .  Results  for  RM  and  Large  Guns 


Quantity 

RM 

Large  Gun 

a 

1.5  x  107m/s2 

1.13  x  106m/s2 

9 .8  cm 

9.9  cm 

a 

vo 

50  volts 

171  volts 

j 

1.92  x  l07A/m 

1.47  x  10?A/m 

Tb 

1.44  x  104  K 

1.83  x  104  K 

t 

51  us 

7.6  ys 

m 

Vm 

0 .88  m/s 

2.3  m/ s 

Vat  xo  = 

=  0) 

114  ys 

419  ys 

tQ  (at  xQ 

=  0.5  m) 

25  ys 

89  ys 

to  ^at  xo 

=  2  m) 

13  ys 

47  ys 

tQ  (at  xQ 

=  4  m) 

9  ys 

33  ys 

t0  (at  xo 

=  5  m) 

8  ys 

29  ys 

The  first  few  entries  in  the  table  were  taken  directly  from  P.ef.  2.  The 

quantities  T,,  t  ,  v  ,  and  tn  were  than  calculated  from  Eqs.  (3.20),  (3.13), 

1  d  m  m  u 

(3,19)  and  (2.4),  respectively.  Several  numerical  calculations  were  undertaken 

8I.R.  McNab ,  "Electromagnetic  Acceleration  by  a  High  Pressure  Plasma ,"  J.  Appl. 

Phys.  Sl_,  2549  (1980). 

9 

S.C.  Rashleigh  and  R.A.  Marshall ,  "Electromagnetic  Acceleration  of  Macropar¬ 
ticles  to  High  Velocities y"  J.  Appl.  Phys.  49_, ,  2540  (1978). 

'  ^ "Laboratory  Demonstration  Electromagnetic  Launcher -EMACK,  "  Commissioning 
Test  Results,  Westinghouse  R&D  Center,  Pittsburgh,  PA. 
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to  make  certain  that  resistive  heating  was  negligible  and  that  Eq.  (3.13) 
adequately  predicted  the  time  to  melting. 

Graphical  results  will  not  be  presented  since  they  are  similar  qualita¬ 
tively  to  those  shown  in  Figs.  2-5.  It  is  interesting,  however,  to  note  that 
for  the  RM  case,  results  similar  to  those  for  the  JB  experiment  are  observed, 
as  can  be  seen  from  Table  II.  Specifically,  the  calculations  predict  rail 
melting  toward  the  breech  end  of  the  gun  (Xq  =  0) ,  but  none  farther  down  the 

gun  tube,  say,  at  Xq  =  0.5m.  For  the  larger  rail  gun,  the  calculated  results 

suggest  melting  throughout  an  entire  5m  gun  tube,  though  the  greatest  amount 
by  far  is  toward  the  breech  end  of  the  gun. 


IV.  SUMMARY  AND  DISCUSSION 

We  have  developed  a  model,  which  when  used  in  conjunction  with  calculated 
or  experimentally  measured  properties  of  the  arc,  can  predict  thermal  damage 
on  the  surfaces  of  the  rails.  Calculations  have  been  undertaken  using  as  in¬ 
put  data  the  experimental  results  of  Jamison  and  Burden,  as  well  as  for  input 
data  obtained  theoretically  in  previous  work.  Results  indicate  substantial 
melting  very  close  to  the  breech  of  the  gun,  with  little  or  none  farther  out 
on  the  gun  tube. 

A  number  of  approximations  are  made  in  order  to  make  the  calculations 
feasible  (see  Sec.  I)  but  in  every  case  the  approximations,  if  anything,  over¬ 
estimate  the  amount  of  wear  on  the  rail  surface.  WTien  using  the  measured  data 
as  input,  the  major  source  of  error  probably  lies  in  the  determination  of  the 
boundary  temperature  of  the  arc  T^.  To  determine  this  temperature  one  must 

rely  on  the  experimentally  measured  muzzle  voltage  V-.  It  has  been  pointed 
8  * 

out  previously,  however,  that  this  measured  voltage  probably  includes  a  sub¬ 
stantial  contribution  from  the  contact  potential  at  the  rail-arc  interfaces. 
Consequently,  the  potential  drop  across  the  arc,  which  should  actually  be  used 
in  Eq.  (3.20),  is  likely  to  be  smaller  than  the  measured  value.  For  instance, 
if  we  assumed  that  only  one-third  of  the  measured  muzzle  voltage  corresponds 

£ 

to  the  potential  across  the  arc,  as  McNab  has  done,  the  model  predicts  for 
the  JB  data  no  melting  even  at  the  point  Xq  =  0.156m.  Again,  the  use  of  too 

large  a  potential  in  the  model  calculations  will  predict  too  large  a  boundary 
temperature  and,  thus,  greater  erosion  than  in  fact  occurs.  The  greatest 
source  of  error  in  using  theoretically  obtained  input  data  is  clearly  in  our 
ignorance  of  accurate  values  of  the  arc  length  %  . 

Uncertainty  in  these  input  data  suggests  the  need  for  additional  work, 
both  theoretical  and  experimental.  On  the  theoretical  side,  two  major  problems 
are  evident.  First,  it  would  be  very  desirable  to  obtain  a  theoretical 
prediction  of  the  arc  length  in  rail-gun  accelerators.  A  study  of  this  prob¬ 
lem  is  difficult  and  would  be  likely  to  involve  a  rather  detailed  understand¬ 
ing  of  how  arcs  are  formed  initially  and  how  they  change  during  the  course  of 
acceleration.  A  second  theoretical  problem  is  to  investigate  the  nature  of 
the  contact  potential  that  no  doubt  exists  at  the  rail -arc  interface.  A  theo¬ 
retical  estimate  of  the  magnitude  of  this  potential  would  clearly  allow  us  to 
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ascertain  what  fraction  of  the  measured  potential  actually  results  from  a  re¬ 
sistive  drop  across  the  plasma. 

From  the  experimental  viewpoint,  more  data  concerning  the  properties  of 
the  arc  at  various  points  along  the  gun  tube  are  much  needed.  Jamison  and 
Burden  have  made  significant  progress  in  the  problem  of  arc  diagnostics,  as 
is  evident  from  the  data  quoted  herein,  but  more  work  remains  to  be  done. 

Some  experimental  analysis  of  rail  erosion  would  also  be  beneficial.  Work  on 
this  problem  has  been  initiated  by  Jamison  and  Burden  using  the  experimental 

nuclear  technique  of  thin-layer  activation.  Some  preliminary  results  have 
been  obtained  but  more  detailed  measurements  are  desirable. 

It  is  hoped  that  as  more  data  concerning  arc  properties  become  available, 
the  model  under  study  here  can  be  used  more  successfully  to  predict  thermal 
damage  to  the  rails.  It  would  also  be  beneficial  in  the  future  to  study  other 
types  of  damage  (say,  mechanical)  both  theoretically  and  experimentally. 


11 
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See ,  e.g.,  A.  Uiiler  and  S.E.  Caldwell,  "The  Fe(p,n)  Co  Reaction  in  Steel 
Wear  Measurement, "  Nucl.  Inst,  and  Methods  12Bj  179  (1976). 


27 


ACKNOWLEDGEMENT 


I  am  grateful  to  J.H.  Batteh  for  several  helpful  discussions,  and  to  K.A. 
Jamison  and  H.S.  Burden  for  providing  me  with  the  data  used  in  the  calculations. 


28 


REFERENCES 


1.  J.D.  Powell  and  J.H.  Batteh,  "Plasma  Dynamics  of  an  Arc-Driven,  Electro¬ 
magnetic,  Projectile  Accelerator,"  J .  Appl .  Phys.  5?,  2717  (1981).  See 
also  "Plasma  Dynamics  of  an  Arc-Driven  Rail  Gun,"~5allistic  Research  Lab¬ 
oratory  Report  No.  ARBRL-TR-02267,  September  1980  (AD  A092345) . 

2.  J.D.  Powell  and  J.H.  Batteh,  "Two-Dimensional  Plasma  Model  for  the  Arc- 
Driven  Rail  Gun,"  J.  Appl.  Phys.  5£,  2242  (1983).  See  also  J.D.  Powell, 
"Two-Dimensional  Model  for  Arc  Dynamics  in  the  Rail  Gun,"  Ballistic 
Research  Laboratory  Report  No.  ARBRL-TR-02423,  October  1982  (AD  A120046) . 

3.  H.G.  Landau,  "Heat  Conduction  in  a  Melting  Solid,"  Quart.  Appl.  Math.  8, 

81  (1950)  .  “ 

4.  P.M.  Morse  and  H.  Feshbach,  Methods  of  Theoretical  Physics  (McGraw-Hill, 
New  York,  1953),  Chap.  7. 

5.  B.  Carnahan,  H.A.  Luther,  and  J.O.  Wilkes,  Applied  Numerical  Methods 
(Wiley,  New  York,  1969),  Chap.  7. 

6.  K.A.  Jamison  and  H.S.  Burden,  "A  Laboratory  Arc-Driven  Rail  Gun,"  Ballis 
tic  Research  Laboratory  Report  No.  ARBRL-TR-02502 ,  June  1983  (AD  A131153) 

7.  K.A.  Jamison  and  H.S.  Burden  (private  communication). 

8.  I.R.  McNab,  "Electromagnetic  Acceleration  by  a  High  Pressure  Plasma,"  J. 
Appl.  Phys.  51_,  2549  (1980). 

9.  S.C.  Rashleigh  and  R.A.  Marshall,  "Electromagnetic  Acceleration  of  Macro¬ 
particles  to  High  Velocities,"  J.  Appl.  Phys.  49^,  2540  (1978). 

10.  "Laboratory  Demonstration  Electromagnetic  Launcher-EMACK,"  Commissioning 
Test  Results,  Westinghouse  R5D  Center,  Pittsburgh,  PA. 

11.  See,  e.g.,  A.  Niiler  and  S.E.  Caldwell,  'The  S^FE(p,n) 56Co  Reaction  in 
Steel  Wear  Measurement,"  Nucl .  Inst,  and  Methods  138,  179  (1976). 


APPENDIX  A 


31 


APPENDIX  A 


In  this  appendix  we  will  be  concerned  with  deriving  an  approximate  ex¬ 
pression  for  the  boundary  temperature  of  the  arc  in  terms  of  easily  measured 
experimental  parameters.  For  purposes  of  calculation,  we  will  neglect  the 
position  dependence  of  the  arc  conductivity  and  radiation  mean  free  path  and 
will  assume  the  condition  w«Jl  .  As  has  been  discussed  before,  neglecting  the 

ci 

position  dependence  in  these  two  quantities  does  not  significantly  affect  the 
results,  whereas  the  second  condition  is  generally  met  in  experiments  of  in¬ 
terest  . 

2 

Under  these  assumptions  the  equation  obeyed  by  the  arc  temperature  is 


(A.l) 


4a  aX  l  2 
s  a 


where  o  corresponds  to  the  arc  conductivity  and  X  to  its  mean  free  path.  In¬ 
tegrating  and  employing  the  obvious  symmetry  condition.  (See  Fig.  1.) 


(A.  2) 


=  0 


we  find 


(A. 3) 


A  second  integration  yields 


(A. 4) 


8a  a  X  Jl 
s  i 


a 


where  is  a  constant. 

Now,  by  assumption,  the  flux  at  y  =  0  satisfies  the  condition 


4 


(A. 5) 


<1  (y-  0)  -  2os  Tb 


However,  if  energy  transport  within  the  arc  occurs  by  radiation  heat  conduc¬ 
tion,  then 


(A. 6) 


q 


3  dy 


Simultaneous  use  of  (A. 4)  -  (A. 6)  then  yields  for  the  boundary  temperature 
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(A. 7) 


.2  1/4 

f _ ^ _ ) 

T  _  v  2J 

b  4a  a  i 


However,  the  potential  drop  across  the  plasma  or  the  "muzzle  voltage”  can  be 
2 

written 


V  =  -2” 
v0  ol 


so  in  terms  of  experimentally  measured  parameters 

jV0  1/4 

Tb  ”  (4o  T"-1 

s  a 


(A.8) 


(A. 9) 
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